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Abstract 



■ Lunar occultation can be used to measure the proper motions of some of the 

o ; 

long time scale microlensing events, t e <; 70 days, now being detected toward the 



Galactic bulge. The long events are difficult to explain within the context of stan- 
dard models of the mass distribution and dynamics of the Galaxy. Han & Gould 
(1995b) have suggested that they may be due to a kinematically cold population 



near the Sun. To resolve the mass, distance, and velocity of individual events and 
so to determine their nature, one must measure parallaxes and proper motions. For 
long events, parallaxes can be often obtained from ground-based measurements, 
but proper motions can only rarely be determined using conventional methods. 
Lunar occultations are therefore key to the understanding of the long events. We 
carry out realistic simulations to estimate the uncertainty of these measurements 
and show that proper motions could be measured for about one long event per 
year. 
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1. Introduction 



The MACHO (Alcock et al. 1995; Bennett et al. 1995) and OGLE (Udalski 
et al. 1994) collaborations have reported 54 candidate lensing events toward the 
Galactic bulge caused by Massive Compact Objects (MACHOs). The current 
estimate of the optical depth obtained by both groups is significantly higher than 
the theoretical estimate (Griest et al. 1991). Various solutions have been proposed 
to explain the observed optical depth excess. One solution assumes a bar-shaped 
Galactic bulge (Kiraga & Pacyhski 1994; Zhao, Spergel, & Rich 1994). However, 
the excess optical depth is not completely solved even with the adoption of a triaxial 
bulge (Han & Gould 1995a). 

Recently, Han & Gould (1995b) argued that the long events, those with time 
scales of t e <: 70 days, cannot be explained within the context of the standard 
models of the mass distribution and stellar dynamics of the Galaxy, leaving the 
origin of these events as a puzzle. Although only three of the reported 54 events 
observed belong to this group, these long events contribute a significant fraction to 
the optical depth, provided that they are truly caused by MACHOs. Hence, their 
potential importance far exceeds their number. Han & Gould (1995b) proposed 
that long events might be caused by a dynamically cold population of dark objects 
such as black holes or neutron stars. Being dynamically cold, such objects would be 
located very close to the Galactic plane, and thus only objects near the Sun could 
contribute to the events detected toward Baade's window located at b = — 3°.9. 
The objects would then have low transverse speeds and so long t e . 

The time scale, t e , is the only observable from current observations, but from 
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t e alone it is difficult to constrain the physical parameters of the long t e events. 
The time scale is related to the physical parameters by 



U = -' r « = c 2 D m ■ (L1) 



Since there are three parameters, and only one measured quantity, the individual 
events are highly degenerate. Here r e is the Einstein ring radius, M is the MACHO 
mass, v is its transverse speed relative to the observer-source line of sight, and D Q \, 
-Dos, and D\ s are the distances between the observer, lens, and source. 

There have been two general approaches to break the degeneracy of the physical 
parameters; parallax and proper motion measurements. If the parallax is measured 
one can obtain the projected speed v: 



v = (1.2) 



The proper motion, /x, is determined by measuring the angular size of the Einstein 
ring, 9 e : 

» = 6 f; 0e = ^. (1.3) 

If both the proper motion and parallax of a MACHO are measured, one can obtain 
the distance to the MACHO by 



A»i = A* (|A* + 1 ) (1-1) 



since the distance to the source stars D os = 8 kpc is known. Similarity M = 
(c 2 /AG)t^vfi and v _1 = v^ 1 + (/xD s) -1 - Therefore, the physical parameters of 



MACHOs would be uniquely determined from the combined information of parallax 
and proper motion. 

For typical short events seen toward the bulge, parallax can be measured only 
from space (Gould 1994, 1995), and proper motions can be measured for only a 
few percent of events where the MACHO transits the face of the star (Gould 1994; 
Nemiroff & Wickramasinghe 1994; Witt 1995). However, parallax measurements 
for long-time-scale events become feasible from the ground while proper motion 
measurement with traditional techniques becomes essentially impossible. Ground 
based parallax measurements are possible because the Earth moves through a sub- 
stantial fraction of the Einstein ring during the event. Indeed Bennett et al. (1995) 
measured a parallax for one of the three long events using only the routine monitor- 
ing data (i.e., no followup photometry). On the other hand, proper motions from 
transits are rare because 9 e <; 1 mas while the angular source size 6* ^ 10 //as, 
even for giants. However, large 9 e opens a possibility of resolving the two images of 
a lensed star and measuring their angular separation, sep . Once sep is measured, 
one can uniquely determine // (see §2). Due to the large 6 e of long events, it may 
one day be possible to measure sep with high-resolution interferometry. Unfortu- 
nately, current interferometry has not achieved high enough resolution to resolve 
the expected very small separations of images; sep ~ 3 mas for long events (see 
§3.2). Instead of interferometry, one can measure sep by using lunar occultation. 
We show in §5.2 one can only measure the angular separation along a certain di- 
rection. However, this uncertainty can often be resolved with information about 
the lens geometry provided by parallax measurements. Even when the information 
from parallax measurement is inadequate, two lunar occultation measurements 
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can resolve the degeneracy. Hence, lunar occultation gives enough information to 
determine the proper motion and so M, v, and D Q \. 

In this paper we analyze the possibility of measuring the angular separation 
of images of a gravitationally-lensed star using the lunar occultation method and 
determine the uncertainty of the measurement by carrying out realistic simula- 
tions. We find that if the observation is carried out for bright, moderate-highly 
magnified Galactic bulge giant stars, one can measure sep with an uncertainty 
^0sep ~ 1 mas. The observational strategy and the method to find the intrinsic 
angular separation from the measured separation normal to the Moon's surface are 
discussed in §5. 

2. Proper Motion From Angular Separation Measurement 

The angular Einstein ring radius and therefore the proper motion of the long 
events can be measured if the angular separation of the two images of a lensed 
star is measured. The locations of the images of the star are the solution to the 
quadratic lens equation given by 

dl _ g sdl _ g2 = 0> ( 2.i) 

where 6\ and #s are the image and source angle measured from the lens location. 
Then the angular separation of the two images is given by 

0sep = \0l+ ~ = (X 2 + 4) 1/2 9 e , (2.2) 

where the dimensionless parameter x is the separation between the source and 
lens measured in units of 9 e : x = 6^/6 e . The value of x is uniquely determined 
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from the light curve because the magnification of the lensed star is related to x 
by A(x) = (x 2 + 2)/x(x 2 + 4) 1 / 2 . Thus, once the angular separation between two 
images of the source star is determined, one can obtain 9 e . 



3. Diffraction Pattern by Lunar Occupation 



3.1. Fresnel Diffraction 

A stellar image produces a Fresnel diffraction pattern when it is occulted by 
the Moon. In this case, the star is idealized as a point source and the Moon's 
disk as a semi-infinite plane. When images of a lensed star are blocked by the 
Moon, the observed diffraction pattern will differ slightly from a perfect point- 
source pattern because the image of the lensed star is composed of two images 
with a small separation. By carefully analyzing the occultation of a lensed star, 
one can use this difference to measure the angular separation which cannot be 
resolved using current telescope technology. Fortunately, the Moon passes through 
the Galactic bulge every month and this permits one to apply the lunar occultation 
method to measure the separation of two images of a star lensed by a significant 
fraction of MACHOs. In addition, since very long events are suspected to be caused 
by MACHOs close to the Sun, the angular separations would be large enough to 
measure sep precisely by the lunar occultation method. 

The diffraction pattern produced by a background source star when it is oc- 
culted by an opaque object (e.g., the Moon) is given by 



where C and S are the Fresnel cosine and sine integrals and go is the intensity 
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(3.1.1) 



without diffraction (e.g., Hecht & Zajac 1979). The two Fresnel integrals are 
defined by 

z z 

C(z) = J cos (^' 2 ) dz 1 , S{z) = J sin (^ /2 ) dz'. (3.1.2) 



Here the dimensionless distance variable z is defined by 

2(doM + duV 1 2 



z = r 



Ac?OMC?M 



(3.1.3) 



where A is the wavelength of the observation, dyi and g?om = -Dos — <^M are the 
distance to the Moon from the Earth and from the source star, and r = dy[9- Since 
doM dM, one can approximate eq. (3.1.3) by z = r(2/ Xdyi) 1 ^ 2 ■ When the image 
of a star is directly on the edge of the Moon's surface, v — r — 0, C(0) =5(0) =0 
and g/go = 0.25. 

3.2. Simulation 

In our simulation the observations are assumed to be carried out as follows. 
During an occultation event, photometry is carried out continuously (i.e., with 
time resolution 1 mas) and read instantaneously. This kind of high speed 
photometry technique has already been developed and the actual instrument (HSP) 
had been installed at the Hubble Space Telescope (HST) although the instrument 
was removed from HST to accommodate the installation of corrective optics. A 
large (^ 4 m) telescope is required for the measurement to compensate for the 
relative faintness of even giant sources in the bulge. For an H = star, a (single 
channel) photometer can detect ~ 9.4 x 10 9 r/AX photons m" 2 s~ 1 /im~ 1 . With 



a band width of AX = 0.3 fim and assuming a detection efficiency rj = 0.5, it 
can detect 1.6 x 10 7 photons ms" 1 using a 4 m telescope with 90 % effective 
surface area. The expected number of photons from a bulge clump giant, typically 
H = 13.2 mag (Tiede, Frogel, & Terndrup 1995), will be ~ 67 ms' 1 assuming 
an extinction of Ajj = 0.25 mag toward the bulge. The sky has a brightness of 
14 mag arcsec~ 2 in the H band producing an expected sky flux of 71 photons ms -1 
in an aperture of 1.5 arcsec in diameter. We assume that the occultation occurs 
when the dimensionless impact parameter x = 0.5, and thus the magnification 
A = 2.18. Then the individual magnifications of the image are 

a x ± Vx 2 + A±x 

A ±= x 2 + _ x 2j X ± = 2 ' ( } 

giving A + = 1.59 and A_ = 0.59 for the primary and secondary images, respec- 
tively. The photon counts for each image are then N Vj \ = 112 photons ms" 1 and 
N Uj 2 = 41 photons ms -1 . The expected number of detected photons from the 
individual images and from the sky for various observational situations are com- 
puted and listed in Table 1. In the table the signal-to- noise ratio is computed by 
S/N = N u /(N Uj i + Ny^) 1 ^ 2 - We assume that other sources of noise, e.g., dark 
current and read-out noise, are negligible. 

The theoretically-expected fringe pattern g/go in the H band, centered at 
1.65 fim, is computed by eq. (3.1.1) and is shown as a function time t in Figure 
1 (a). In the figure, time is measured from the moment when the first (closer to 
the Moon) image just crosses the Moon's limb and negative time implies that the 
image is behind the Moon. The angular distance is related to time scale by 

6 = utcosi), (3.2.2) 



where u is the orbital angular speed of the Moon and ip is the angle of the lunar 
limb relative to the direction of the Moon's motion. The total number of photons 
in the signal is normalized into unity when t — oo. Due to the phase shift of the 
Fresnel integrals with time, there are beat patterns. In the simulation, we assume 
the angular separation between the two images of lensed stars is sep = 3 mas, 
equivalent to an Einstein radius 9 e ~ 1.5 mas, which would be typical for a disk 
lens located at ~ 2 kpc with mass of M = 0.7 M@. Since the long events are 
probably caused by relatively massive MACHOs, the expected separation could be 
greater. 

Up until now, our computation has been based on a monochromatic wave 
observation. However, the if -band filter has a finite bandpass. For the correction 
of band width in our analysis, we average the flux weighted by the filter function. 
The filter function in the H band is well-approximated by the tophat function 
^filter: 



where uq = (0.3 //m) -1 . The resultant fringe pattern h(9), after being averaged 
by the filter function, is shown in Figure 1 (b). The small-scale fluctuation in flux 
is smeared out, and thus the beat patterns, which might have been useful for the 
measurement of (f> se p, disappear. 

We make the simulation more realistic by including the beam pattern of the 
telescope. At any time, one edge of the telescope mirror will see an image displaced 
by a small angle compared to the image seen at the center of the mirror. The 
resulting image is then the combination of all images seen by different parts of the 




ujq, 1-5 /im < A < 1.8 /dm 
0, otherwise, 



(3.2.3) 
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mirror. The beam pattern decreases with increasing angular separation 9 P from 
the center of the mirror because the relative area of the mirror in a strip located 
at 9 p decreases as 9 P increases. For a circular mirror /beam = [l — (#p/#tel) 2 ] , 
where 9 te \ = a/2dyi = 1 mas for an assumed aperture a = 4 m. Then the expected 
intensity h{9) is obtained by convolving the intensity h{9) with the beam pattern 

/beam($p) : 

1(0) = J /beam(0p " 9)h(9 p )d9 p . (3.2.4) 

The resulting fringe pattern after the beam pattern convolution and being averaged 
by filter function is shown in Figure 1 (c). 

4. Uncertainty of Measurement and Observational Strategy 

The determination of the angular separation of the two images of a lensed 
star is made by fitting the observed fringe pattern to a set of light curves with 
different S ep- This fitting process suffers from another uncertainty in addition to 
the one from the angular separation of the images. This uncertainty comes from 
the ambiguity of a reference point, which is the time of first image occultation, (i.e., 
9 = 0). A misalignment of the reference point would result in misinterpretation 
of sep . Therefore, we include the misalignment quantified by the shift in the 
reference point, Ato, as a free parameter in our analysis. 

The contours of equal uncertainties measured by x 2 i n the parameter space 
of Ato an d 4>sep are shown in Figure 2 (a) for the first lensing event described in 
Table 1. The contour levels are drawn at 1 a, 2 a, and 3 a levels from the point of 
minimum % 2 marked by "x" . In the contour map there appear two minima which 
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are located symmetrically about sep = 0. These two minima appear because 
of the ambiguity of order of occultation: the fringe pattern when the primary 
image approaches the Moon first looks similar to the pattern when the other image 
approaches first. If the two images have exactly the same intensity, there will be no 
difference in the resultant fringe patterns. The uncertainty is, A(p sep ~ 0.55 mas 
at the 1 a level. This would be a significant detection although not as precise as 
one might like. 

However, the uncertainty decreases significantly with increasing S/N ratio. 
There are two ways to increase the S/N ratio; observing bright stars or highly- 
magnified events. For illustration, we compute the uncertainties for events under 
exactly the same conditions as the previous case except for a higher magnification 
A = 3.0 [Fig. 2 (b)] and a 0.5 mag brighter source star [Fig. 2 (c)], which are the 
second and third cases in Table 1, respectively. For both cases the uncertainty is 
Z\0 sep ~ 0.35 mas at the 1 a level. 

5. Practical Considerations 

5.1. Lunar Topography 

For an actual observation, one is required to consider some miscellaneous fac- 
tors that make sep deviate somewhat from what we have assumed in the compu- 
tation. The first deviation arises because the Moon's limb is not a perfect straight 
line but has a small curvature. If the angular separation between images is big 
enough (e.g., binary stars), the observed pattern would slightly deviate from the 
Fresnel diffraction pattern, which assumes an infinite straight-line surface. How- 
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ever, for the scale of & 5.5 m, which is equivalent to sep ^ 3 mas projected on the 
Moon's surface, the Moon's limb is a nearly perfect straight line. Indeed, typical 
errors arising from lunar surface curvature are negligible even for the measurement 
of a few tens of arcsec (e.g., Richichi et al. 1994). The other type of deviation 
which does affect on the measurement of sep is due to topographical features on 
the lunar surface. If the images are occulted at different parts of a geological struc- 
ture on the Moon, e.g., moutains and valleys, the angle ip in eq. (3.2.2) deviates 
from what one assumes based on a smooth circular lunar disk. Assuming that the 
deviation due to large scale strucures of order km or larger can be corrected using 
a detailed lunar map, one still has a problem due to smaller-scale structures such 
as rocks and cliffs. Fortunately, lunar craters and other topographical features are 
generally not very rugged, and even mountains are very smooth (Abell, Morrison, 
& Wolff 1993). 

5.2. Orientation of Images 

We have assumed up until now that the orientation of the two images is perpen- 
dicular to the approaching limb of the Moon. However, the orientation in general 
is random and thus what one measures is the component of sep normal to the 
Moon's limb not the intrinsic separation of interest. If one knows exactly how the 
source star crosses the Einstein ring, it is simple geometry to deduce sep from the 
measured angular separation (pi, where the subscript "i" is explained below. How- 
ever, resolving the orientations of the images is not a trivial problem and requires 
additional information. 

Information about the lens geometry can be obtained from (ground-based) 
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measurement of the parallax. At a minimum, parallax measurements determine 
the component of the projected velocity parallel to the ecliptic and the magnitude 
of the component normal to the ecliptic. In general, this still leaves a four-fold 
degeneracy: two-fold for the sign of the normal component and two-fold for the 
sense of source motion (clockwise or counterclockwise) relative to observer-lens line 
of sight. In some cases, especially, if an event is observed away from the ecliptic 
and has a long time scale, this four-fold degeneracy is completely broken by the 
parallax measurement itself. Indeed, this is the case for the long-event parallax 
measured by Bennett et al. (1995). In these cases the geometry is unambiguous 
and 9 e can be determined directly from sep . However, it is difficult to break this 
degeneracy when the event is close to the ecliptic. Since the Moon never gets 
farther than ~ 5° from the ecliptic, it is also important to consider the degenerate 
case. 

The four-fold degeneracy can be broken by carrying out occultation observa- 
tions twice so that the approaching angle, 9 a, of the Moon's limb is different for 
each observation. This can be done either by observing the occultation at two dif- 
ferent locations on the Earth or by observing an occultation from the same location 
a month later. The latter case is possible for long events since the declination of 
Moon's orbit typically changes by ~ 1/3 of its diameter per month at fixed right 
ascension. The four-fold degeneracy of lens geometries is illustrated in Figure 4: 
transverse velocities pointing below and above the ecilptic are noted by "case I" 
and "case II". Subnotations, "(a)" and "(b)", describe the clockwise and counter- 
clockwise source crossings. For illustration, we assume that the Moon's limb is 
approaching the images with 9 a = 45° ("/" sense) and 9 a = —45° ("\" sense). 
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For each observation shown in the figure the shaded surface indicates the orien- 
tation of the surface of the Moon and the shaded arrow indicates the direction of 
motion. The angular separation that is measured when 9a = 45° and 6>a = —45° 
is marked by 4>i an d 02- For each case the ratios between measured angular sepa- 
rations determined at two different locations (or times), 0i/02, are different from 
one another. Therefore, one can obtain the full lens geometry by ruling out the 
other three cases from a comparison of the observed ratio 0i/02 with the expected 
ratios. 

5.3. Event Rate 

One can measure proper motions using lunar occultations for about one long 
time-scale lensing event per year. The occulation can be observed in a strip of sky 
along the path of the Moon with a width larger than the angular size of the Moon 
by observing the occultation at higher (or lower) Earth latitudes, 0. The effective 
occulting cross-section is ~ 1°.5 assuming the observations could be arranged at 
places within = ±30°. Then the area of sky toward the bulge, with width ±5° 
around the Galactic plane, swept by the effective cross-section is 15 deg . If the ob- 
servation could be carried out at more extreme Earth latitudes, the event rate could 
be increased. The MACHO group has detected 13 giant events out of a total of 44 
events during a bulge season by covering 12 deg of sky. Three of 44 are long events. 
The expected occultation event rate is then 3 x (13/44) x (15/12) ~ 1 event/yr. 
This is substantially higher than the number of proper motions of long lensing 
events that can be obtained using other methods. Another advantage of the oc- 
cultation method is that for the long events the measurement could be repeated 
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for better determination of proper motion. Special attention should be paid to 
observing those fields which the Moon will occult sometime during a given bulge 
season. For some regions close to the plane, optical observations are impossible 
due to extinction. However, events can still be detected using H or K band obser- 
vations. Indeed, the long events may well be concentrated near the plane (see Fig. 
1 from Bennett et al. 1995). Because the long events last many months, it would 
be sufficient to make these infrared observations once per week compared to once 
per day for standard optical observations. 

Acknowledgement: We would like to thank M. Everett and G. Newsom for 
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